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' This chapter provides a review of microscopic theories of pairing in nuclear sys- 

, terns and neutron stars. Special attention is given to the mean-field BCS theory 

and its extensions to include effects of polarization of the medium and retardation 
of the interactions. Superfluidity in nuclear systems that exhibit isospin asymme- 
try is studied. We further address the crossover from the weak-coupling BCS 
^ ' description to the strong-coupling BEC limit in dilute nuclear systems. Finally, 

^ \ within the observational context of rotational anomalies of pulsars, we discuss 

P5 . models of the vortex state in superfluid neutron stars and of the mutual friction 

between superfluid and normal components, along with the possibility of type-I 
superconductivity of the proton subsystem. 



1. Introduction 

Neutron stars represent one of the densest concentrations of matter in our universe. 
These compact stellar objects are born in the gravitational collapse of luminous stars 
with masses exceeding the Chandrasekhar mass limit. The observational phenomena 
characteristic of neutron stars, such as the pulsed radio emission, thermal X-ray 
radiation from their surfaces, and gravity waves emitted in isolation or from binaries, 
provide information on their structure, composition, and dynamics. The properties 
of superdense matter are fundamental to our understanding of nature at small 
distances characteristic of nuclear forces and of the underlying theory of strong 
interactions - QCD. In fact, neutron stars (NS) provide a unique setting in which 
all of the known forces ~ strong, electroweak, and gravitational - play essential roles 
in determining observable properties. 

Except for the very early stages of their evolution, neutron stars are extremely 
cold, highly degenerate objects. Their interior temperatures are typically a few 
hundreds of keV, far below the characteristic Fermi energies of the constituent 
fermions, which run to tens or hundreds of MeV. Although very repulsive at short 
distances, the strong interaction between nuclcons is sufficiently attractive to induce 
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a pairing phase transition to a superfluid state of neutron-star matter. As will be 
discussed in this chapter, the existence of superfluid components within neutron 
stars has far-reaching consequences for their observational manifestations. 

Historically, the first observational evidence for superfluidity in neutron-star in- 
teriors was provided by the timing of radio emissions of pulsars, the first class of 
neutron stars, discovered in 1967 by Jocelyn Bell. The pulsed emissions, with a 
typical periodicity of seconds or less, are locked to the rotation period of the star. 
Although pulsars are nearly prefect clocks, their periods increase gradually over 
time, corresponding to a secular loss of rotational energy. Significantly, some pulsars 
are found to exhibit deviations from this impressive regularity. The pulsar timing 
anomalies divide roughly into three types, (i) Glitches or Macrojumps. These are 
distinguished by abrupt increases in the rotation and spin-down rates of pulsars by 
amounts Aft/n ~ 10"^- IQ-^ and A(l/n ~ IQ-^. After a glitch, An/n and An/n 
slowly relax toward their pre-glitch values, on a tim e scale of order weeks to years, 
in some cases with permanent hysteresis effects.^^Such behavior is attributed to a 
component within the star that is only weakly coupled to the rigidly rotating nor- 
mal component responsible for the emission of pulsed radiation - an interpretation 
supported by fits of the measured rotation under different modeling assumptions. 

(ii) Timing Noise or Microjumps. These represent irregular, stochastic deviations 
in the spin and spin-down rates that are superimposed on the near-perfect periodic 
rotation of the star. The origin of microjumps remains unclear, but they could be 
evidence of stochastic coupling between the superfluid and normal components.El 

(iii) Long-Term Periodic Variabilities. Observed in the timing of few pulsars, most 
notably PSR B1828-11, these deviations strongly constrain theories of superfluid 
friction inside NS, if their periodicities are interpreted in terms of NS precession.l^ 

X-ray observations from orbiting spacecraft, which yield estimates of surface 
temperature for a half-dozen or so yo ung neutron stars, further reinforce the picture 
of NS with superfluid content! I ^ l ^ l < ' At the stellar ages involved, neutrino emission 
from the dense interior dominates thermal evolution, with nucleonic superfluidity 
acting to suppress the main emission mechanisms. The existing measurements of 
surface temperatures indicate that superfluid hadronic components must be present 
in some NS, since otherwise they would cool to temperatures below the empirical 
estimates on very short time scales. Finally, since NS with their huge gravitational 
fields are expected to be major sources of gravitational wave radiation, it is believed 
that observation of gravitational waves from oscillating neutron stars can provide 
further information on the state of matter in their interiors. In particular, the eigen- 
frequencies and damping rates of gravity waves may carry imprints of dissipation 
processes in the superfiuid phases^S'9 Vj* 

The existence of neutron-star superfluidity, first envisioned by Migdal m 
1959, is broadly consistent with microscopic theories of nucleonic matter in NS. 
Shortly after the advent of the Bardeen-Cooper-Schrieffer (BCS) theory in 1957, 
BCS pairing of nucleons in nuclei and infinite nuclear matter was suggested and 
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Fig. 1. Left panel. Dependence of experimental scattering phase shifts in ^S\, ^Pi, ^D2, and 
partial waves on laboratory energy. Right panel. Dependence of critical temperatures of superfluid 
phase transitions in attractive channels on chemical potential. The corresponding densities are 
indicated by arrows. 



studied! ' With the discovery of puls ars, the imphcations of nucleonic pairing 
for neutron-star properties were explored}^ a.nd viable micr oscopic calculations of 
pairing gaps began to appear soon thereafter! -*^ ^'^I'*^^ ! -*^^ ! -*^^ ! '*^'^ ! 

Partial-wave analysis of the nucleon-nucleon {NN) scattering data yields infor- 
mation on the dominant pairing channel in nuclear- and neutron-matter problems in 
a given range of density (see Fig.^. At high densities, corresponding to laboratory 
energies above 250 MeV, the most attractive pairing channel is the tensor-coupled 
'^P2~^F2 channelj'^O ^I^whenever isospin symmetry is even slightly broken. This con- 
dition holds inside neutron stars, with the partial densities of neutron and proton 
fluids differing quite significantly, except in special meson-condensed phases where 
the nucleonic matter is isospin-symmetric. In such a case the phase-s hift ana lysis 
predicts that the most attractive pairing interaction is in the •^I?2 wave.E3^ 

At low density, isospin-symmetric nuclear matter exhibits pairing due to the 
attractive interaction in the ^Si-'^Di partial wave, a tensor component of the force 
again being responsible for the coupling of the S and D waves. '^^ '^'^ This 
interaction channel is distinguished by its ability to support a two-body bound 
state in free space - the deuteron. However, under the highly isospin-asymmetric 
condition typical in neutron stars, neutron-proton pair condensation is quenched by 
the large discrepancy between the neutron and proton Fermi momenta (see Scc.|21l. 

This review is devoted to several aspects of nucleonic superfluidity in neutron 
stars that are of major current interest. Section |2] provides an overview of the many- 
body theories of pairing in neutron stars, with a special focus on the Green's function 
description of pairing, and the effects of self-energies and vertex corrections. The 
key ideas of the correlated-basis (or "CBF" ) approach to superfluid states of Fermi 
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systems are also presented. Section 13 is concerned with the possibiUty of pairing 
between particles lying on different Fermi surfaces, in particular, between protons 
and neutrons in asymmetric nuclear matter. Here we determine the critical value of 
isospin asymmetry at which a transition from '^Si-'^Di to ^5*0 pairing occurs. (This 
threshold value is small compared to p — n asymmetries typical of neutron-star 
matter). More generally, we consider several competing phases that could exist in 
asymmetric mixtures of fermion species. There exist systems of this kind, notably 
dilute, ultracold atomic gases and baryonic matter as described by QCD at high 
density, where "asymmetric pairing" is enforced. (In the case of dilute atomic gases, 
such conditions can be "tuned in" by external fields, while in the QCD problem, 
the conditions of charge neutrality and f3 equilibrium together with the heaviness 
of the strange quark lead naturally to asymmetric pairing among light quarks.) 
In Sec. 21 we turn to the phenomenon of crossover from BCS superconductivity to 
Bose-Einstein condensation, as it occurs in fermionic systems that support a two- 
body bound state in free space. Section [S] is devoted to the physics of superfluids 
at the "mesoscopic" scale, with discussions of flux quantization, neutron vorticity, 
the electrodynamics of superconducting protons, and their implications for mod- 
eling the dynamics of rotational anomalies in pulsar timing. Our conclu sions a nd 
related perspectives are summarized in Sec. |S| Two other recent reviewJ^SEHI of- 
fer complementary information and perspectives on nuclear pairing and nucleonic 
superfluidity. 

2. Many-body theories of pairing 
2.1. Propagators 

This subsection outlines the Green's functions method for the treatment of super- 
fluid systems. The original formulations of this approach are due to Gor'kov and 
Nambu, who employed thermodynamic Green's functions!^ Herein we consider 
the real-time, flnite-temperature formalism, which is suited to studies of both equi- 
librium and non-equilibrium systems. Our discussion is restricted to equilibrium 
systems, and we shall work with the retarded components of the full Green's func- 
tion of the non-equilibrium theory. Superfluid systems are described in terms of 
2x2 matrices of propagators (known as Nambu-Gor'kov matrix propagators) that 
are deflned as 



where ipa{x) are the baryon field operators, x is the space-time coordinate, the 
indices a and /3 stand for the internal (discrete) degrees of freedom, and T and T 
denote time-ordering and inverse time-ordering of operators, respectively. The 2x2 
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matrix Green's function satisfies the Schwinger-Dyson equation 



/ , I a. u, .i, 

7,(5 

where the free-propagator matrices Gapi^^^') ^tre diagonal in the Nambu-Gor'kov 
space. The matrix structure of the self-energy flapiXjX') is identical to that of the 
propagators: the on-diagonal elements are and p), and the off-diagonal ele- 
ments are A(p) and A'^{p). Fourier transforming Eq. ^ with respect to the relative 
coordinate x — x' , one obtains the Dyson equation in the momentum representa- 
tion. In this representation, the components of the Nambu-Gor'kov matrix obey the 
coupled Dyson equations 

Gc^pip) = Goaisip) + Goajip) [^^s{p)Gsf3{p) + A^s{p)f1 spip)] , (3) 

F^pip) = Goaj{-p) [Als{p)Gsp{p) + ^,s{~p)F^sp{p)] ■ (4) 

Here Gap{p) and Goa/^ip) are the full and free normal propagators, F''ap{p) and 
Fap{p) are the anomalous propagators, and Y^a/iijp) and Aq^(p) are the normal 
and anomalous self-energies. The Greek subscripts are the spin/isospin indices, and 
summation over repeated indices is understood. For systems with time-reversal sym- 
metry, it is sufficient to solve Eqs. Q and since this symmetry implies that 
Aafj{p) = [A^^(p)]*. It is instructive to rewrite Eqs. in terms of auxiliary 

Green's functions 

G^p{p) = GoaMp) + G^^{p)i:,s{p)Go5p{p) (5) 

describing the unpaired state. The solution of this equation is G^p{p) = 6ai3[uj — 
e{p)]~^, where e{p) = ep + and €p is the free single-particle spectrum. (N.B. 
Assuming that the forces conserve spin and isospin, the self-energy E(p) is diagonal 
in spin and isospin spaces). Combining Eqs. ©,(01), and ISJ, we derive an alternative 
but equivalent form of the Schwinger-Dyson equations, namely 

GM = G^^ip) [d,p + A^s{p)F^s0ip)] , Flf,{p) = G^^{-p)A\sip)Gspip) , 
which can be solved to obtain 



[u: - EA{p)f - Es{p)^ - A^p) 



pUp) - 7 — ^...2"'::.o .o. . • (7) 



Here we have made the substitution A(p)At(p) = — A^(p) and defined symmetric 
and antisymmetric parts of the single-particle spectrum in the normal state, Eg^j^ = 
[sip) ± e{—p)] /2. For isotropic systems, the self-energy is invariant under reflections 
in space (i.e. the self-energy is even under p — > —p); furthermore, for systems 
that are time-reversal invariant, the self-energies are even under the transformation 
Lu —uj. Hence the antisymmetric piece of the spectrum Ea must be absent when 
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both conditions are met. The poles of the propagators © determine the excitation 
spectrum of the superfLuid system, given by 

Lo± = EAip)±VEsip)^ + ^Hp)- (8) 

Here one sees that there is a finite energy cost ~ 2A for creating an excitation from 
the ground state of the system when = 0, a property that leads to the existence 
of superflow or supercurrcnt in paired fermionic systems. 

The solutions © and Q are completely general, all functions being dependent 
on the three-momentum and the energy. Superfluid systems are often treated in 
the quasiparticlc approximation, in which the self-energies are approximated by 
their on-mass-shell counterparts. The rationale behind such an approach is that the 
nuclear system in its ground state can then be described in terms of Fermi-liquid 
theory, if one neglects pair correlations. Switching on the pair correlations precipi- 
tates a rearrangement of the Fermi-surface, but it is assumed that the quasiparticle 
concept remains intact. 

Within this framework, the wave-function renormalization is defined by ex- 
panding the normal self-energy as S(u;) = ^{£p) + di^T,{uj)\^^i.^{uj — Sp), where 
Ep = tp -|-ReI](ep) is the on-mass-shell single-particle spectrum in the normal state 
that solves Eq. Now, it is seen that the propagators © and Q retain their 
form if they are renormalized as 

G{p)=Z{p)G{u + i5,p), F^p) = Z{p)F\cj + tS,p), K\p) = Z {pf ^\p) , 

(9) 

where Z{p) = [l — di^Ti{ijj)\^=^J\ is the wave-function renormalization and the 
tilde identifies renormalized quantities. An additional feature of the renormalized 
propagators is that the quasiparticle spectrum e(p) is now constrained to the mass 
shell. (N.B. For time-local interactions the gap function is energy- independent, so 
there is no need to expand A(ix') around its on-shell value. We shall return to the 
off-shellness of the self-energies in Subsec. [T^ \ 

Renormalization of propagators within the quasiparticle picture suggests that 
the probability of finding an excitation with given momentum p is strongly peaked at 
the value £p. The wave- function renormalization takes into account corrections that 
are linear in the departure from this value. From the computational point of view, 
such corrections require a knowledge of the off-shell normal self-energies. The depen- 
dence of the imaginary part of the self-energy (quasiparticle damping) on frequency 
follows from Fermi-liquid theory, being given by 2ImI](a;) — a 1 -I- (cj/27rT)^ , 
where a is a density-dependent constant. The real part of the self-energy, ReS([jj), 
can be computed from ImS(Li;) via the Kramers-Kronig dispersion relation only if 
the latter function is known for all frequencies uj e [— oo, oo\. Accordingly, the fore- 
going result from Fermi-liquid theory should be supplemented by a model of the 
high-energy tail of the quasiparticle damping. 

The momentum dependence of propagators can be approximated by introducing 
an effective quasiparticle mass. For nonrelativistic particles, expansion of the normal 
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self-energy around the Fermi momentum leads to 

1 -1 

m 

■ (10) 



f -s PF f ^ * ™ 

£[P) = —[P- PF) ~ fJ- , 

TO* TO 



TTL 

1 + — 5pReI](p)|p=p^ 

PF 



where fi* = —e{pF) + fJ- — Rel](pi?). A closed system of equations determining the 
properties of the superfluid system is obtained by specifying the self-energies in 
terms of the propagators and interactions, as discussed below. 

2.2. Mean- field BCS theory 

The BCS-type theory of superconductivity as applied to nuclear systems is predi- 
cated on a mean-field approximation to the anomalous self-energy. Specifically, the 
anomalous self-energy is expressed through the four-point vertex function T(j),p') 
in the form 

A(P) = -2 1 ^nP.P') l^F{p')f{u:') , (11) 

where f{u!) — [1-f exp(/9cj)]~^ and (3 is the inverse temperature. We observe that the 
gap is energy-independent when the interactions are local in time, corresponding 
to no retardation, as in the case where the effective interaction T{p,p') is replaced 
by the bare interaction V{p,p'). In this case, carrying out th e renorm alization 
according to Eq. and integrating over the energy, we arrive aJ^S^ESI 

A(p) = Z{p) I J^Vip,p')Zip')^^ [/(C.+ ) - /(c._)] , (12) 



Aip) = Z{p) I -l^V{p,p')Z{p').r^^ [/(^+) - /(^_)] , (13) 



where uj± = ±^ + IS? [cf. Eq. 0]. Further progress requires partial- wave decom- 
position of the interaction in Eq. (|12|l . To avoid excessive notation, our consideration 
focuses on a single, uncoupled channel, for which we obtain a one-dimensional gap 
equation 

dpP^.rf f^yf nMpl 
{2ttY U}+[P) 
where V{p,p') is the interaction in the given partial wave. The gap equation is 
supplemented by the equation for the density of the system, 

"'-'oj Ifi'-'^w/i-) = i /^^w _ i:_ (i + 1) /(-) . (") 

which determines the chemical potential in a self-consistent manner. Here g is the 
isospin degeneracy factor. Eqs. (|11() - (|14|1 define mean-field BCS theory. In subse- 
quent discussions based on this theory, the tilde notation will be suppressed. 

Assuming that the interaction V{p^p') and the single-particle spectrum Sp in 
the unpaired state are known, Eqs. ((T^ and l(TH) form a closed system for the gap 
and chemical potential. Within this formulation, the presence of a hard core in the 
interaction causes no overt problems, for one can readily show that at the critical 
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temperature, Eq. transforms into an integral equation that sums the particle- 
particle ladder series to all orders. In the special case in which the interaction is 
momentum- independent, Eq. (|13l) leads to the familiar weak-coupling formula 



where the density of states is specified by j^(pf) ~ m*pFZ^{pF)/2TT^ (for one di- 
rection of isospin). The weak-coupling formula is often used to estimate the mag- 
nitude of the gap. However, because realistic pairing interactions are momentum- 
dependent, the density dependence of the gap predicted by Eq. may deviate 
markedly from that given by the weak-coupling formula H15(l . Moreover, as argued 
in Ref . KMI if the pairing interaction acquires a strong momentum dependence due to 
a short-range repulsive core, the weak-coupling formula may well produce a mean- 
ingless or useless estimate of the gap, since its derivation requires that V{p,p') take 
a negative value on the Fermi surface. 

The normal-state self-energy is written as 



where the subscript A indicates antisymmetrization of the final states and the con- 
tribution cx ImT is omitted for simplicity. For nuclear systems, the amplitude T is 
often approximated by the scattering T-matrix, which sums up the ladder diagrams 
and is generally nonlocal. An alternative is to replace the T-matrix by an effective 
time-local interaction that is fitted to properties of finite nuclei (e.g. a Skyrme or 
Gogny force), in which case Eq. IjlGI) reduces to a mean- field Hartree-Fock approx- 
imation. However, one must be aware that the normal-state spectrum itself must 
depend on the anomalous self-energy A(p). The common replacement of G{p) by 

(p) when computing the normal-state spectrum is an approximation (sometimes 
called the "decoupling approximation"), which is justified when the pairing effects 
can be viewed as a perturbation to the normal state, but must be made with care. 

Fig.|21shows the pairing gap in neutron matter and symmetrical nuclear mat- 
ter for the high-precision phasc-shift-equivalent Nijmegen potential and the eff ecti ve 
Gogny DSl force, for different approximations to the single-particle spectrum.!^ In 
the case of smooth effective forces such as the the Gogny interaction, a Hartree- 
Fock approximation to the normal self-energy is suitable and was adopted; for the 
(realistic) Nijmegen potential, the T-matrix was calculated in Brueckner theory. In 
both cases, the full momentum-dependent self-energies Rel](p) were used in the gap 
equation. The momentum renormalization yields an effective mass m* /m less than 
unity, thus reducing the density of states ^{pf) on the Fermi surface and hence also 
reducing the size of the gap. The wave-function renormalization factor Zij)) is also 
less than one, leading to an additional suppre ssion o f the gap. However, the the 
magnitude of this effect is yet to be established P^^^ 




(15) 




(16) 
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Fig. 2. Singlet 5-wave (^So) pairing gap in neutron matter and nuclear matter versus Fermi mo- 
mentum. The heavy and light lines correspond respectively to the phase-shift-equivalent Nijmegen 
interaction and the Gogny effective interaction DSl. Solid lines and dashed-dotted lines label 
results for pure neutron matter using the free-space single-particle spectrum and implementing 
single-particle renormalization, respectively, while the dashed line refers to symmetrical nuclear 
matter with single-particle renormalization. 

2.3. Polarization effects 

An improvement upon the mean-field BCS approximation to fermion pairing is 
achieved in theories that take into account the modifications of the pairing interac- 
tion due to the background medium. In diagrammatic language, the class of modi- 
fications known as "polarization effects" or "screening" arise from the particle-hole 
bubble diagrams, ideally summed to all orders starting from the bare interaction 
as the driving term. Consider the following integral equation describing the four- 
fermion scattering process pi + p2 ^ Ps + Pi- 

r(p, p',q)^U{p,p',q)-lJ^ U{p, p", q)G^ {p" + q/2)G^{p" - q/2)T{p",p', q) , 

(17) 

where q = pi~ P2 is the momentum transfer, p = pi +p3, and p' ~ P2+PA- Eq- l|17|l 
sums the particle-hole diagrams to all orders. To avoid double summation in the 
gap equation, the driving term U{p^p',q) must be devoid of blocks that contain 
particle-particle ladders. This driving interaction depends in general on the spin 
and isospin and can be decomposed as 



Uq = fq + 9q{a- ■ ct') + [/^ + .g^(cr • cr')] (r • t') , 



(18) 
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where cr and r are the vectors formed from the Pauh matrices in the spin and 
isospin spaces. We assume here that the interaction block U depends only on the 
momentum transfer. For illustrative purposes, the tensor part of the interaction 
and the spin-orbit terms are ignored. (However, see Subsec. 12.41 where the tensor 
component is included by means of pion exchange.) Solution of the integral equation 
(|17|l then takes the form 

Go 



F„ 



l + A(g)Fq l + A(g)G. 



-{a 



F' 



l + K{q)F' l + A(g)G: 



where Fq = v{pF)fq, Gq = v{pF)9q, F' 



^{PF)fq, and G'q 



{t-t'), 

= ^{PF)9'q, while 



K{q) = v(p,. 



(2^ 



G^(y' + (7/2)G^(p"-g/2) 



(19) 



(20) 



is the (dimensionless) Lindhard function (or polarization tensor). We are tacitly 
assuming that the system is in a state characterized by a well-defined Fermi sphere. 
Then, if the momenta of both particles lie on the Fermi surface, the momentum 
transfer is related to the scattering angle via q — 2pi?sin0/2, and the parameters 
F^ F' , G, and G' can be expanded in spherical harmonics with respect to the 
scattering angle, according to 



Fiq)\ 
G{q)) 



E 



Pi (cos ( 



(21) 



and similarly for F'{q) and G'{q). The Landau parameters Fi, Gi, Fl, and G[ depend 
only on the density. The isospin degeneracy of neutron matter, reflected in x-r' — 1, 
implies that the number of independent Landau parameters for each q or I reduces 
from four to two, defined by F" = F + F' and G" — G + G' . Commonly, only the 
lowest-order harmonics in the expansion (|21() are needed. For a singlet pairing state, 
in which the total spin of the pair is 5 = and cr ■ cr' = —3, the pairing interaction 
is given by 



iy{pF)T, = FS 



A{q)F^ 



3G^ 



A(g)Gg 
1 + A(9)GS 



(22) 



1 + A{q)F^ 

In general, the polarization tensor A{q) is complex- valued. However, in the limit of 
zero energy transfer (at fixed momentum), it is real and becomes simply 



A(g) 



T 



4pl 



In 



2pi 



2pi 



q 



(23) 



in the zero-temperature limit. Eq. H22|l contains two distinct contributions: the 
direct part generated by the terms 1 inside the square brackets, and the remaining, 
induced part that accounts for density and spin-density fluctuations. If the Landau 
parameters are known - either by inferring them from experiment or by computing 
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Fig. 3. Singlet 5-wave (^5o) pairing gaps in neutron matter versus Fermi momentum, as obtained 
from microscopic calculations that attempt to account for medium-modification of the pairing 
interaction and self-energies. The curves are labeled as: A - Wambach et al. [42], B - Chen et al. 
[48], C - Chen et al. [49], D - Schwenk et al. [52], E - Schulze et al. [51], F - Fabrocini et al. [50]. 



them within an ab initio many-body scheme - the effect of polarization can be 
assessed by defining an averaged interaction 

1 r'^PF 

r = -^/ dqqT{q). (24) 

The effect of density fluctuations cx is to enhance the attraction in the pair- 
ing interactiorP^ and therefore increase the gap, while the spin-density fluctua- 
tions (X Gg tend to reduce the attraction and decrease the gap^H At densities 
typical of the inner crust of a neutron star, the values of microscopically derived 
Landau parameters imply th at t he suppression of pairing via spin-density fluctu- 
ations is the dominant effect!^ The Landau parameters of neutron matter and 
symmetrical nuclear matter have been studied extensively within complex many- 
body scheme d I I I whose description is beyond the scope of this chapter. It 
should be noted, however, that the matrix elements of the pairing interactions de- 
rived in some of these schemes, including the B abu-Br own approach"^" and its 
extensions such as polarization-potential theory, | 41|42 | j^g^yg i^gen used in conjunc- 
tion with the weak-coupling approximation (|15|l . which - as indicated above - is 
generally inadequate in the nuclear/neutron-matter context. 

The initial microscopic calculations of the effects of medium polarization on 
pairing were carried out within an a.lternative many-bod y approach, the method 
of Correlated Basis Functions (CBF; ' ■ ^ 1 I ^ ■ ' 1 ^ 1 i l ^v l ^V I {^q be described in Sub- 
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sec. I2.5f) . At a qualitative level, the findings of the CBF studies of Chen et a 
are consistent with much of the later work based on Green's functions and Fermi- 
liquid theories. Until recently, there was broad agreement that the screening reduces 
the singlet 5- wave gap by factor of 3 or so (c.f. Fig. [21). However, the density profiles 
of the calculated gaps differ considerably. This is illustrated in Fig.O which presents 
a composite plot of theoretical predictions for the dependence of the ^5*0 pairing gap 
at the Fermi surface, A{kp), upon the Fermi wave number kp = pp/h = (Stt^p)^/"^. 
The six curves in the plot correspond to various microsc opic approaches that in- 
clude a screeni ng c orrection. The results of Wambach et al.,E21 Schulze et al.,l^and 
Schwenk et al.l^ are based on microscopic treatments rooted in Landau/Fermi- 
liquid theorVjWith polarization-type diag rams summed to all orders. The results of 
Chen et al.ESEni and Fabrocini et al.l^ were obtained within two different imple- 
mentations of CBF theory. 

Further assessment of the status of quantitative microscopic evaluation of the 
singlet- S" gap is deferred until Subsec. [2.51 where the elements of CBF approaches 
to the pairing problem are reviewed. We may already remark, however, that explicit 
comparison of the pairing matrices constructed in the different theories could help 
to eliminate discrepancies introduced by use of the weak-coupling approximation in 
some of the theoretical treatments. Another important consideration is consistent 
inclusion of medium effects on both the pairing interaction and the self-energies. 

2.4. Non-adiabatic superconductivity 

Since mesons propagate in nuclear matter at finite speed, the interactions among 
nucleons are necessarily retarded in character. As a consequence, the self-energies 
(and in particular the gap function) must depend on energy or frequency. Within the 
meson-exchange picture of nuclear interactions, the lightest mesons - pious - should 
be the main source of nonlocality in time. This suggests that it may be fruitful to 
consider a pairing model in which the interactions are modeled in terms of pion 
exchanges, plus contact terms that can be approximated by Landau parameters. 
Thus one assumes an interaction structure 

Vnn = - — {(T-V){T-cl}) + U{q) (25) 

in which cp is the pseudoscalar isovector pion field satisfying the Klein-Gordon 
equation, is the pion-nucleon coupling constant, and m^r is the pion mass. Here 
(f) is the pseudoscalar isovector pion field satisfying the Klein-Gordon equation, is 
the pion-nucleon coupling constant, and TOtt is the pion mass. The operator structure 
of the term U (q) is like that of Eq. H18|l , but with constants differing numerically 
since the tensor one-pion exchange is treated separately. 

For static pions, the one-pion-exchange two-nucleon interaction in momentum 
space is given by 

K(.) = - J, ^ kr • + SMn)] rr . r. , (26) 
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Fig. 4. Top panel. Baryon Fock-exchange self-energies for normal (left graph) and anomalous 
(right graph) sectors. The solid lines correspond to fermions; and the wavy lines, to pions. The 
blobs symbolize RPA-renormalized vertices; and the dots, bare pion-nucleon vertices. Bottom panel. 
Normal and anomalous Hartree diagrams. The dots stand for contact Landau interactions; the 
straight lines arc shown for clarity. 

where q is the momentum transfer, n = q/q, and Si2{n) is the tensor operator. This 
interaction is known to reproduce the low-energy phase shifts and, to a large extent, 
the deuteron properties.''^ Below, however, the pairing correlations are evaluated 
from the diagrams that contain dynamical pions, with full account of the frequency 
dependence of the pion propagators. The static results can be recovered, and a rela- 
tion to the phase shifts established, only in the limit ^ in the pion propagator. 
The dynamics at intermediate and short range is dominated in turn by the cor- 
related two-pion, p-meson, and heavier-meson exchanges. Short-range correlations 
are crucial for a realistic description of low-energy phenomena, being necessary to 
achieve nuclear saturation. Moreover, the response functions calculated from one- 
pion exchange alone would already precipitate a pion-condensation instability in 
nuclear matter at an unrealistically low density. 

The interaction (|25|l leads t o a time- nonlocal formulation of nuclear supercon- 
ductivity in neutron-star matter.l^ The Dyson-Schwinger equation (O provides the 
starting point. Since the formulation will incorporate the full frequency dependence 
of the self-energies, it is useful (and also conventional to define the wave- function 
renormalization differently than in Subsec. 12.21 Thus we set Z{p) = 1 — 
the retarded self-energy being decomposed into components even (S) and odd (A) 
in uj, i.e. = "^sip) + The single-particle energy is then renormalized as 

Es — + T,s{Es,p). Accordingly, the propagators now take the forms 



G{p) 



wZ{p) + Es{p) 



(27) 



F{p) 



{uo + irifZipf - Es{pY - A(p)2 ' 
A(p) 



(28) 



(w + iry)2Z(p)2 - EsipY - A(p)2 ' 



where AA^^ = 



— A2. The self-energies of the theory are shown in Fig. 2| The ana- 
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lytical counterparts of the Fock self-energies are 



E'^-"^(u;,p) = -2 p^^Eo{q)lmG{s,p-q)C{L.,s,q)Eiq), (29) 



{2n) 

A^-^a;,p) = -2 p^E,{q)lTnF {e , p - q)C {u , e , q)^{q) , (30) 



where 



rooi / 

C{i^,e,q)= -^B{u;',q) 



f{e)+g{uj') , l-f{e) + g{uj') 



IT] e + LU' — LU — M] 



(31) 



Here B{q) is the pion spectral function, while S°((j) and S(g) are the bare and 
renormalized pion-neutron vertices. One remarkable feature of Eqs. (|29|I - 13H) is 
that the energy and momentum dependence of the self-energies is determined by 
the dynamical features of the meson (here pion) field. Another salient feature is that 
the normal and anomalous sectors are coupled, in contrast to the BCS case, where 
the unpaired single-particle energy is unaffected by the pairing. The most impor- 
tant contribution to the pion spectral function comes from the coupling to virtual 
particle-hole states, which are described by the (retarded) particle-hole polarization 
tensor Il{uj,p). Specifically, one finds 

^ ^ -2Imn^(g) 

[cj2 _ q2 _ _ Ren«(q)]2 + [Inin«(q)]2 ' ^ ' 

The spectral function of pions in neutron matter is illustrated in Fig. [S] for the 
momentum transfer range < g < 2kp, where kp = 0.55 fm"-'^. It is seen that 
(i) the spectral function has a substantial weight for finite energy transfer, the max- 
imum being determined by the pion dispersion relation uP = q^ + + Ren(a), q) 
and (ii) the spectral function is substantially broadened due to the excitations of 
particle-hole pairs, which are treated in the random-phase approximation (RPA).^3 
In addition to the Fock self-energies we need to include the Hartree contribution, 
which reduces to 

A^^CS(^ p)^_2(Fo"-3GS) j 0^lmF{c. + c.',p + p'). (33) 

Solutions of the self-consistent equations H30|l and H31II are shown in Fig. at 
zero temperature and for densities specified by the indicated Fermi momenta. At 
small energy transfers, the imaginary components of the gap and wave- function 
renormalization vanish, and one recovers the BCS limit. For finite energy transfers 
these functions develop substantial structure that refiects the features of the pion 
spectral function [the driving term in the kernel of integral equations for A{q) and 
Z{q)]. Note that the actual value of the gap on the mass shell does depend on 
the detailed structure of these functions far from the mass shell. However, it is 
possible to renormalize the pion spectral function such that the high-energy tails 
are eliminated while the on-mass-shell physics is unchanged. 
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Fig. 5. Pion spectral function in neutron matter as a function of energy and momentum transfer. 
The density corresponds to kp = 0.55 fm~^. 

Theories that expHcitly include the hght mesons - pions or kaons - in the com- 
putational scheme have the advantage that they embody the precursor phenomena 
associated with the softening of the pion (kaon) modes close to the threshold for con- 
densation. In the c ase of P-wave pairing, an enhancement of the pairing correlation 
has been predicted.!^ 



2.5. The method of correlated basis functions ( CBF) 

The variational formulation of BCS theory is based on the trial wave function (BCS 
state) 

*Bcs = n [(1 - V)'^' + ^y'V^tV'lpx] |0) ' (34) 

p 

where the real function hp, which gives the occupation probability of the pair state 
[p Ti ~P i)j is subject to variation. For this trial state one may compute the anoma- 
lous density 

Xp = (^BCsl^T^-Pil^BCs) = hl/'{l - hpfl^ , (35) 

whose Fourier image specifies the spatial structure of a Cooper pair. The standard 
coupled BCS equations for the energy gap, the fermionic density, and the quasipar- 
ticle energy at zero temperature, namely 



A{p) 



(27r)3 



(36) 
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Fig. 6. Top panel. Frequency dependence of real (heavy lines) and imaginary (light lines) parts 
of the gap function, Ai(a;) and A2{u!), respectively, for kp = 0.4 (solid), 0.5 (long-dashed), 0.55 
(dashed), and 0.6 (dashed-dotted lines) fm~^. The on-shell values of the pairing gap are 0.1, 0.7, 
1.4, and 3.7 MeV for kp = 0.4, 0.5, 0.55, and 0.6 fm~^. Bottom panel. Frequency dependence of 
real (left panel) and imaginary (right panel) parts of the wave-function renormalization, Z\{ll)) 
and Z2(t<j). Labeling is the same as in the top panel. 



(2^)3^ 



p^dp 
(2^ 



2E{p) 



(37) 



E{p) = d[e{p)~^iY + ^^{p). 



(38) 



are generated naturally upon (i) evaluating the expectation value of the grand- 
canonical Hamiltonian for the BCS state (|34|) . and (ii) performing a variational 
minimization of this functional with respect to hp, under the constraint that the 
expectation value of the particle-number density coincides with the prescribed den- 
sity. 

Treatment of strongly interacting fermionic systems (including nuclear prob- 
lems) within a variational framework calls for improved supcrfluid trial states. The 
systems of interest are characterized by a bare two-body interaction containing a 
strong inner repulsive core along with longer-range attractive components. To obtain 
a reasonable energy expectation value, the trial function must adequately describe 
the short-range geometric correlation s induced by the repulsion, which inhibits the 
close approach of a pair of particles ! ^ v l ^^ l ^ v l ^O I xhe simplest choice involves the 
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Jastrow correlation factor 

N 

Fj = n /('-y ) ' 1™ /(O ^ , lim /(r) ^ 1 , (39) 

i<j 

which is suited to efhcient description of state-independent two-body correlations, 
especially the short-range repulsive effects. As a bonus, with proper optimization of 
the two-body function /(r), the Jastrow factor can also incorporate effects of virtual 
phonon excitations and indeed can reproduce the correct asymptotic behavior of 
long-range correlations. 

A substantially improved trial superfluid state of definite particle number may 
be formed by applying the Jastrow operator l|39|l to an A'^-particle projection of the 
BCS trial state, expressed in the configuration-space representation as 

•^■1303 = ;^-4{</>(ri,r2),0(r2,r3)...0(r7v-i,rw)}. (40) 

Here, 0(ri,rj) is the antisymmetrized Fourier image of Xp given by Eq. H35|l . times 
a spin function, A is an antisymmetrization operator acting on particle pairs in dif- 
ferent 4> factors, and is a normalization constant. In spirit, this is the approach 
adopted in the early work of Yang and Clark,^ although in practice they de- 
termined the one-body, two-body, etc. density matrices required for their cluster- 
expansion treatment of the short-range correlations from the BCS state (|34|l . rather 
than from its A^-particle projection (|4(J|I . 

Following up on the work of Refs. El I16L a formal variational theorjEHl of the 
superfluid ground state of uniform, infinite nucleonic systems was developed for a 
trial correlated BCS state constructed in Fock space, 

1*^) -EE F^''^^(^^){'^iV\QCS) , (41) 

where F'^'^^ is an unspecified correlation operator meeting certain minimal condi- 
tions and {l^m'^)} is a complete set of Fermi-gas Slater determinants, both referred 
to the A^-particle Hilbert space. Thus, the correlated normal states {F^^'^\<^m'^)} 
are superposed with the same amplitudes as the model states 1$^^) have in the 
corresponding grand-canonical representation of the original BCS state (|34|l . Rep- 
etition of steps (i) and (ii) above for this correlated superfluid Ansatz yields a 
theory having the same structure as ordinary BCS theory, when the "decoupling 
approximation" is applied. In this approximation, only one Cooper pair at a time 
is considered, while treating the background as normal. Formally, the expectation 
value of the grand-canonical Hamiltonian is expanded in terms of the deviations of 
the Bogolyubov amplitudes Up = {1 — hpY^'^ and Vp = h]/^ about their normal- 
state values, retaining deviant terms at most of flrst order in Vp — 6{p) and second 
order in UpVp, where 9{p) is the Fermi step. 

Within this framework, the gap equation and density constraint maintain the 
same mean-field forms as obtained for the bare BCS state, except for the attach- 
ment of renormalization factors z^^ to the gap function A(p) when it appears in 
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quasiparticle energy denominators. The presence of correlations introduced by the 
operator i^(^) is otherwise reflected only in the replacement of the pairing matrix 
elements V{p,p') and single-particle energies s{p) derived from the bare interaction 
based on the BCS trial state 134(1 and variational steps (i)-(ii), by effective pairing 
matrix elements V{p,p') and correlation-dressed single-particle energies ^(p) built 
from combinations of diagonal and off-diagonal matrix elements of the Hamiltonian 
and unit operators in the correlated normal bases {F*^^^ j*!)^'')}. When the state- 
independent Jastrow choice Fj of Eq. H39|) is assumed for the correlation operator 
F, the dressed quantities V{p,p') and £^{p) can be evaluated by Fermi hypernetted- 
chain (FHNC) methods developed in Ref. El with results for neutron matter and 
liquid "^He reported in Ref. 

An important advance in the CBF approach to pairing was made in Ref. 1471 
where the variational description was extended to create a correlated-basis pertur- 
bation theory for the exact superfluid ground state. Again imposing the decoupling 
approximation, a sequence of approximations to the grand-canonical energy may be 
defined, each preserving under variation the standard form of the gap equation, but 
with successive improvements on the effective pairing matrix elements and dressed 
self-energies. (A convenient modification of the trial correlated ground state (|41|l 
was made by inserting the normalization factor (^J^-* jF^^^^i^'^^ j'l'^^'')^^/^ inside 
the summations. This eliminates the renormalization factors mentioned above.) 
Making the Jastrow choice for the correlation operator F, the leading perturbative 
corrections to the variational results for V{p,p') and ^{p) were generated, repre- 
sented in diagrammatic form, and evaluated. These corrections include the leading 
contribution from medium polarization within the CBF framework. Here we should 
point out that the terms in the CBF perturbation expansions of the various quan- 
tities are not easily interpreted in terms of conventional Goldstone or Feynman 
diagrams, although they may have similar appearance. A given order in the CBF 
expansion, including the "zeroth-order" variational term, will contain pieces be- 
longing to any number of perturbative orders in the conventional sense. In general 
there will be terms accounting for the nonorthogonality of the basis, terms that cor- 
rect the average-propagator approximation inherent (for example) in the Jastrow 
description of F-correlations, terms that correct for non-optimality of the chosen 
F-correlations, etc. 

In microscopic studies of nucleonic systems, it is generally imperative to include 
the effects of state-dependent correlations arising from realistic NN interactions 
which contain, separately in each spin and isospin channel, contributions of central, 
tensor, and spin-orbit character. In principle, the CBF perturbation expansions 
provide for systematic correction of the Jastrow Ansatz for the correlation operator 
F, so as to incorporate these state-dependent effects. However, it is clearly preferable 
to take account of state dependence already in the choice of F, thus reducing the 
need to correct the variational treatment with CBF perturbation theory. A suitably 
general correlation Ansatz^ within the class containing only two-body correlation 
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factors, is given by 



N 



F 



^ij ) Oct 



(42) 



i<j a—1 



where f{ij) contains terms for the same operators Oa{ij) as are present in the as- 
sumed reahstic NN interaction (e.g., the Argonne uig model"^^), or an adequate 
subset of them. Profound difficulties arise in the implementation of this choice, due 
to non-commutativity among the Oa{ij) operators. The analog of FHNC resum- 
mation being still beyond our reach for such state-dependent correlations, existing 
calculations proceed with st raightf orward cluster or power-series expansions, per- 
haps with vertex corrections! I I It is important to appreciate that the extended 
Jastrow form H42() of the i^-operator is equipped to include the lion's share of the po- 
larization corrections (just as the simple state-independent Jastrow form is capable 
of capturing the major effects of density-density fluctuations). 

Chen et al.ESl applied CBF pairing theory as developed in Ref. H7ito superfluid 
neutron matter in the phase, assuming state-independent Jastrow correlations 
and taking account of the leading CBF perturbation corrections to the variational 
treatment. The polarization and other corrections produced a very substantial sup- 
pression from the variational estimate of the gap A(fci?), the peak value being re- 
duced by a factor ^ 4 and situated at lower density. This treatment was updated in 
Ref. Major improvements were made in the choice of the variational two-body 
correlation functions. Significantly larger gap values were obtained, but again the 
perturbative corrections were estimated to suppress the peak value by a factor ~ 4 
and shift its location to a lower density. Quantitatively, the results of the later of 
the two perturbative CBF calculations are the more reliable. 

Shortly after the work of Krotscheck and Clarkj^Hl an independent approach to 
CBF description of pairing based on the trial superfluid state H41|l was launched 
by Fantoni.l^ With immediate specialization of the correlation operator F to the 
state- independent Jastrow form H39|l . it proved feasible to extend the diagrammatic 
techniques of standard Fermi hypernetted-chain theorj^^ and thereby enable 
practical evaluation of the one-body density matrix and radial distribution function 
associated with the correlated BCS state (I41II . Derivation of the corresponding gap 
equation and density constraint was achieved without resorting to the decoupling 
approximation. 

Quite recently,!^ Fantoni's CBF approach has been generalized - insofar as 
practicable - to include state-dependent correlations of the form H42I) . The results for 
the ^Sq neutron gap lie distinctly higher than the results of earlier work designed to 
include nontrivial medium effects on pairing (see Fig.lJJ. Thus, they conflict with the 
general consensus that these effects lead, on balance, to a strong suppression of the 
gap value. Concurrent estimateJ^ of the gap based on an auxiliary-field diffusion 
Monte Carlo (AFDMC) calculation lie even higher than the new CBF estimates (by 
roughly 0.5 MeV, with a peak value of more than 2.5 MeV ai kp 0.6 fm^^). A 
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recent numerical study of ^5*0 pairing in neutron matter within the self-consistent 
Green's function (SCGF) methoa^''' gives results in essential agreement with the 
AFDMC estimates. One might conclude from this agreement that the medium- 
polarization effects, arising from the exchange of spin-density fluctuations and other 
virtual processesj^are less important than previously imagined. On the other hand, 
the SCGF calculation, by construction, neglects such collective correlations of longer 
range, while the AFDMC stochastic estimates, obtained for relatively small samples 
of neutrons, might also fall short in their inclusion of these effects. At any rate, the 
latest computational results continue to highlight the extreme sensitivity of the 
^So pairing gap to the assumptions made in pursuing its evaluation by microscopic 
methods. Thequantitative situation for pairing in spin-triplet T = 1 states is even 
less clear.ESES 



3. Pairing in asymmetric nuclear systems 

The isospin asymmetries characteristic of neutron-star cores, with proton fractions 
~ 5%, are too large to permit isospin-singlet (neutron-proton) pairing. A pos- 
sible exception involves Bose-Einstein condensation of kaons at densities several 
times nuclear saturation density, in which case the matter is approximately isospin- 
symmetric. In high-density is ospin-s ymmetric nuclear matter, neutron-proton pairs 
form in the partial wave.^^ZEHl jjowever, once the isospin symmetry is slightly 
broken, pairing is suppressed and isospin-triplet neutron-neutron and proton- 
proton pairs are formed. Due t o their large partial density, neutrons pair in the 
^P2-^F2 tensor -coupled channelpSEl while the less abundant protons pair in the 
state.EEni 

We have seen in Subsec. (12.1(1 that for fermionic systems which are invariant un- 
der reversal of time and reflections of space, the quasiparticle spectrum is symmetric 
under p —p, and consequently the antisymmetric piece Ea in Eq. vanishes. 
Depending on the system, these symmetries could be broken either by the presence 
of external gauge fields or due to intrinsic properties such as the mass difference in 
mixtures of gases. At any rate, we now focus on systems having Ea ^ 0, and the 
pairing in question is between fermions that lie on different Fermi surfaces. We shall 
call such systems asymmetric superconductors(hereafter ASC). 

Initial studies of ASC where carried out in the early sixties when, shortly af- 
ter development of the BCS theory of superconductivity, metallic supe rconductors 
with paramagnetic impurities were studied experimentally! I I I ' Since colli- 
sions with impurities can flip the spins of electrons, an imbalance between spin-up 
and spin-down electron populations is created. This effect can be mimicked by in- 
troducing an average, effective magnetic field that lifts the electron-spin degeneracy 
due to its interaction with the electron magnetic moment. The novel aspect of the 
studies of ASC (apart from the new context) is the realization that a correct in- 
terpretation of the results requires a self-consistent solution of the gap and density 
equations, even in the weak-coupling limit where the changes in the value of the 
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Fig. 7. Left panel. Temperature dependence of pairing gap for density asymmetries a = 0.0 
(solid), 0.05 (daslied-dotted), 0.07 (dashed), and 0.1 (dotted). Right panel. Temperature depen- 
dence of free energy. Labeling of asymmetries is as in the left panel. 



chemical potential due to pairing are small. 

To avoid undue complications in describing ASC, we will operate within the 
framework of conventional BCS theory, in the sense that effects of wave-function 
renormalization and medium polarization are neglected. We shall also suppress the 
additional complication of the ^^i-^Di tensor coupling. This aspect is not essential 
for the present discussion; see Refs. 1241 1251 1261 for the relevant details. 

Thus, the equations underlying the theory of ASC are taken to be and H14|l . 
the spectrum being given by Eq. (|SJ) with 7^ 0. If the spatial symmetries are 
unbroken, then Ea = Sjjt = /ip)/2, where /i„ and /ip are the neutron and proton 
chemical potentials. In general, Eqs. (|13ll and H14f) must be solved self-consistently. 
Consider first the procedure in which 

Eq. is solved by parametrizing the asymmetry in terms of the difference 
in the chemical potentials, and the densities of the speci es are computed after the 
gap equation is solved. Such an analysis predict J^^^^^^ a double- valued character 
of the gap as a function oi Sji. On the first branch, the gap has a constant value 
A((5/x) = A(0) over the asymmetry range < i5/Lt < A(0) and vanishes beyond the 
point Sji — A(0). The second branch exists in the range A(0)/2 < (S/i < A(0), 
with the gap increasing from zero at the lower limit to A(0) at the upper limit. 
Only the portion 5^ < A{0)/^/2 of the upper branch is stable, i.e., it is only in 
this range of asymmetries that the superconducti ng s tate lowers the grand thermo- 
dynamic potential from that of the normal state.l^ Thus, the dependence of the 
superconducting state on the shift in the Fermi surfaces is characterize d by a con- 
stant value of the gap, which vanishes at the Chandrasekhar-Clogston^^^ limit 
Sfii = A(0)/V2. 

A different picture emerges from an alternative treatment of the problem in 
which particle- number conservation is incorporated explicitly by solving Eqs. 1)13(1 
and (|14l) self-consistently.*^^ -^ These studies find a single- valued gap as a function 
of the isospin asymmetry a = {pn — Pp)/{Pn + Pp)- Minimizing the free energy of 
an asymmetric superconductor at fixed density and temperature leads to stable 



February 9, 2008 0:39 WSPC/Trim Size: 9.75in x 6.5in for Review Volume Sedrakian 



22 A. Sedrakian and J. W. Clark 

solutions over the entire reRion of density asymmetries where non-trivial solutions 
of the gap equation exist M'^"'^ This can be seen in Fig. |3 where the temperature 
and asymmetry dependence of the pairing gap and the free-energy of a homoge- 
neous asymmetric superconductor are shown. In particular, we see that for a fixed 
temperature, the gap and the free energy are single-valued functions of the density 
asymmetry a in the particle-number-conserving scheme - in contrast to what is 
found in the non-conserving scheme, where double-valued solutions appear. 

At large asymmetries, the dependence of the gap on the temperature shows 
a "re-entrance" phenomenon. As the temperature is increased from low values at 
which the asymmetry is too large to sustain a gap, a critical temperature is reached 
at which pairing correlations take hold. (For example, this is seen for the a = 0.1 
case in Fig.[7|). This behavior can be attributed to the increase of phase-space over- 
lap between the quasiparticles that pair, due to the thermal smearing of the Fermi 
surfaces. Further increase of temperature suppresses the pairing gap at a higher 
critical temperature due to thermal excitation of the system, in much the same 
way as in the symmetric superconductors. Clearly, in this scenario the pairing gap 
has a maximum at some intermediate temperature. The values of the two critical 
temperatures are controlled by different mechanisms. The superconductivity (or 
superfluidity) is destroyed with decreasing temperature at the lower critical tem- 
perature when the smearing of the Fermi surfaces becomes insufficient to maintain 
the required phase-space coherence. The upper critical temperature is the analog 
of the BCS critical temperature and corresponds to a transition to (re-entrance of) 
the normal state because of thermal excitation. 

Another as pect of the asymmetric superconducting state is the gapless nature of 
the excitationsl^^^l^ One may draw an analogy to the non-ideal Bose gas, for which 
only a fraction of the particles are in the zero-momentum ground state at temper- 
atures below the critical value for Bose condensation. The dynamical properties 
of gapless superconductors, such as response to electroweak probes and transport, 
depend on the ratio ( = A/ Sfi in an essential way: for C > 1, the response of the 
system is similar to that of an ordinary superconductor; in the opposite limit C < !> 
the system's behavior is essentially non-superconducting (see e.g. Ref . 17^ . These 
features are easily understood by examining the excitation spectrum in both limits. 



3.1. Phases with broken space symmetries 

We now turn to a special class of ASC characterized by broken global symme- 
tries - translational, rotational, or both. In 1964, Larkin and Ovchinnikov—' and, 
independently, Fulde and Ferrell^ (LOFF) discovered that the superconducting 
state can sustain asymmetries beyond the Chandrasekhar-Clogston limit if electrons 
pair with nonzero center-of-mass (hereafter CM) momentum. The weak-coupling 
result for the critical shift in the Fermi surfaces for onset of the LOFF phase is 
6fi2 — 0.755 A(0) [> 5fii = 0.707 A(0)] . Since the condensate wave function depends 
on the CM momentum of the pair, its Fourier transform will vary in configuration 
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Fig. 8. Left panel. Dependence of pairing gap A in the LOFF phase on density asymmetry a and 
total momentum P of the condensate, relative to pairing gap Aqo in the limit P = = a. Right 
panel. Dependence of free energy of the LOFF phase on the same quantities as in the left panel. 

space, giving rise to a lattice structure. The Fulde-Ferrell state is predicated on a 
plane-wave form A(r) = Aexp(— zP • r) for the gap function. Larkin and Ovchin- 
nikov considered a number of lattice types and concluded that the body-centered- 
cubic (bbc) lattice is the most stable configuration near the critical temperature. 
Recent studies of this problem in the vicinity of the critical temperature employ- 
ing Ginzburg-Landau theory show that the face-centered-cubic (fee) structure is 
favored.l^ 

For illustrative purposes, let us consider the Fulde-Ferrell state, in which case 
the quasiparticle spectrum is given by 



where the upper sign corresponds to neutrons and the lower sign to protons. The 
spectrum (I43II is obtained by applying the following transformation to Eq. JH)): 
Es ^ Es + (P2 -t- p^)/2m and Ea Sfi ± p ■ P. Onset of the LOFF phase 
entails a positive increase in the quasiparticle kinetic energy cx Q^, which disfavors 
the Fulde-Ferrell state relative to the BCS state. However, the anisotropic term 
oc P-p, which can be interpreted as a dipole deformation of the isotropic spectrum, 
modifies the phase-space overlap of the fermions and promotes pairing. The LOFF 
phase becomes stable when the increase in the kinetic energy required to move the 
condensate is smaller than the reduction in potential energy made possible by the 
increase in the phase-space overlap. The magnitude of the total momentum serves as 
a variational parameter for minimization of the ground-state energy of the system. 
The pairing gap and the free energy of ASC with finite momentum are shown in 
Fig. IS] It is assumed that the gap function depends paramet ric a lly on the magnitude 
of the CM momentum, but is independent of its direction.C^ For such an Ansatz 
the anisotropy of the spectrum appears only in the Fermi functions in the kernels 
of Eqs. H13I) and (|14|l and is averaged through the phase-space integration. It is 




(43) 
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seen in Fig. |S| that an ASC-LOFF state arises for arbitrary finite momentum of the 
condensate below some critical value. For large enough asymmetries the minimum 
of the free energy moves from P = to intermediate values of P, i.e., the ground 
state of the system corresponds to a condensate with nonzero CM momentum of 
Cooper pairs. Note that for the near-critical range of asymmetries, the condensate 
exists only in the LOFF state and its dependence on the total momentum exhibits 
the re-entrance behavior found in the temperature dependence of the homogeneous 
ASC. The order of the phase transition from the LOFF to the normal state is a 
complex issue that depends on the preferred lattice structure, among other things 
(see Ref. 1771 and work cited therein). 

Due to the Pauli exclusion principle, a noninteracting fermionic gas fills an 
isotropic Fermi sphere; similarly, if there are two types of noninteracting fermions, 
each species fills an isotropic Fermi sphere. Consider now a strongly interacting 
system that is a Fermi liquid rather than a Fermi gas. According to Fermi-liquid 
theory, the states of the interacting system are reached by switching the interac- 
tion on adiabatically. Driven by this process, the noninteracting gas evolves into a 
strongly-interacting liquid, in which the dressed single-particle degrees of freedom 
- the quasiparticles - once again fill a spherical shell isotropically. However, this 
simple Fermi-liquid picture may not hold in two-component (or multi-component) 
fermionic systems in which the fermions of differing species interact via strong pair- 
ing forces. Indeed, there can exist a stable superconducting phase that sustains 
ellipsoidal deformations of the Fermi-s urfaces , a phase hereafter referred to as de- 
formed Fermi-surface 

superconductivitjI^SIZni (DFS). 

The quadrupole deformations of the Fermi surfaces are described by expand- 
ing the quas ipartic le spectrum in spherical harmonics and keeping the I — 2 
contributions) ' v l ( ? ! 



where i^±{q) is the spectrum of the homogeneous ASC and the coefficients £2 de- 
scribe the deformations of the Fermi surfaces that break the rotational 0(3) sym- 
metry down to 0(2). The 0(2) symmetry axis is chosen spontaneously. Thus, the 
quasiparticle spectrum of the DFS phase is obtained from the spectrum of homo- 
geneous ASC by the transformations 



We observe that the leading harmonic term responsible for deformation of either 
Fermi surface is that for 1 = 2, not 1 = 1, since the latter corresponds to a translation 
of one Fermi sphere relative to the other, without deformation. The deformations 
are deemed to be stable if they lower the free energy of the system relative to its 
value in the undcformed state. 

The deformation parameters e2± are determined by minimization of the free 
energy of the system, as was done for the Cooper-pair momentum parameter P 
in the case of broken translational invariance (LOFF states). Three-dimensional 




(44) 



Es^ Es + (£2,+ + e2,-)/2M , Ea^Ea + (£2,+ - e2,-)/2A* • 



(45) 
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Fig. 9. Left panel. Dependence of pairing gap in the DFS phase on density asymmetry and total 
momentum of the condensate. Right panel. Dependences of free energy of the DFS phase for the 
same input parameters as in Fig. 5. 
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Fig. 10. Projection of the Fermi surfaces on a plane parallel to the axis of symmetry breaking. 
The concentric circles correspond to the two populations of spin(isospin)-up and spin(isospin)- 
down fermions in a spherically symmetric state (<5e = 0), while the deformed symbols correspond 
to the state with relative deformation St = 0.64. The spin{isospin)-density asymmetry is a = 0.35. 

plots of the dependence of the pairing gap and free energy of the DFS phase on 
asymmetry and the relative deformation Se = (e2,+ ^ e2,-)/2/i are provided in 
Fig. |5| Fig. ^1 shows a typical deformed Fermi-surface configuration that lowers 
the expected ground-state energy below that of the non-deformed state. At a = 0, 
the critical deformation for which pairing ceases is the same for prolate and oblate 
deformations. At finite a and in the positive range of Se, the maximum value of the 
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gap is attained for constant 5e; at negative Se the maximum increases as a function 
of deformation and saturates for (Se ~ 1. The re-entrance phenomenon sets in for 
large asymmetries as Se is increased from zero to finite values. (N.B. The essential 
difference between LOFF and DFS phases is that in the latter, the translational 
symmetry of the superconductor remains unbroken.) 

To complete our discussion of pairing states in nonrelativistic asymmetric su- 
perconductors, we briefly mention some of the alternatives to the LOFF and DFS 
phases. One possibility is that the system prefers a phase separation of the super- 
conducting and normal phases in real space, such that the superconducting phase 
contains particles with matching chem ical potentials, i.e. is symmetric, while the 
normal phase remains asymmetric 

Equal-spin (-isospin, -flavor) pairing is another option, if the interaction between 
like-spin particles is attractive^ Since the separation of the Fermi surfaces does 
not affect spin-1 pairing on each Fermi surface, an asymmetric superconductor may 
evolve into a spin-1 superconducting state (rather than a non-superconducting state) 
as the asymmetry is increased. Therefore spin-1 pairing becomes the limiting case 
for very large asymmetries. If the single-particle states defining the different Fermi 
surfaces are characterized by spin (as is the case in the metallic superconductors), 
the pairing interaction in a spin-1 state should be P-wave and the transition is 
from S-wave to P-wave pairing. If the fermions are characterized by one or more 
additional discrete quantum numbers (say isospin as well as spin), the transition 
may occur between different S'-wave phases (e.g. from isospin-singlet to isospin- 
triplet in the case of nuclear matter). The possibilities become especially rich in 
dense quark matter. 



4. Crossover from BCS pairing to Bose-Einstein condensation 

A crossover from BCS superconductivity to Bose-Einstein condensation (BEG) is 
exhibited in fermionic systems with attractive interactions under sufficient decrease 
of the density and/or sufficient increase of the interaction strength. The transition 
from large overlapping Cooper pairs to tightly bound non-overlapping bosons can 
be described entirely within the ordinary BCS theory, if the effects of fluctuations 
are ignored (mean- field approximation). Early studies of this type of transition were 
carried out in the contexts of ordinary superconductors, excitonic superconduc- 
tivity in semiconductors,!^ and, at finite temperature, an attractive fermion gas.l^ 
Although the BCS and BEC limits are physically quite different, the transition 
between them is found to be smooth within ordinary BCS theory. 

Several authors have considered the BCS-to-BEC transition in the nuclear con- 
text. In isospin-symmctric nuclear matter, neutron-proton (np) pairing undergoes 
a smooth transition leading from an assembly of np Cooper pairs at higher den- 
sities to a gas of Bo se-condensed d euterons as the nucleon density is reduced to 
extremely low values! I < | oo | o^ | b)U | 'Yhis transition may be relevant - and could 
then yield valuable information on np correlations ~ in low-density nuclear sys- 
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Fig. 11. Dependence of pairing gap (upper panel) and chemical potential (lower panel) on tem- 
perature for fixed values of the ratio / = no/n, where n is the baryon density and jiq = 0.16 fm~^ 
is the saturation density of symmetrical nuclear matter. Values of the diluteness parameter na'^ 
assume a scattering length a = 5.4 fm. 



terns (especially the nuclear surface), in expanding nuclear matter from heavy- 
ion collisions, and in supernova matter. The underlying equations of the theory 
are (113(1 and ((14(1 with Ea = 0; we shall address the effects of asymmetry at a 
later point. The top panel of Fig. II II shows the dependence of the gap function on 
temperature for several densities given in terms of the ratio / = no/n, where 
no = 0.16 fm~'^ is the saturation density of symmetrical nuclear matter. The bottom 
panel shows the associated chemical potentials /i computed self-consistently from 
Eq. II14I) . The low- and high-temperature asymptotics of the gap function are well 
described by the BCS relations A(r ^ 0) = A(0) - [27raA(0)r]i/2 exp(-A(0)/T) 
and A(r Tea) = 3.06 /3[Te2(Te2 - T)]i/2, respectively, where T^a is the criti- 
cal temperature of the phase transition. However, the BCS weak-coupling values 
a — 1 — f3 must be replaced with a ~ 0.2 and (3 ~ 0.9. As a consequence, the 
ratio of the zero-temperature gap to the critical temperature deviates from the fa- 
miliar BCS result A(0)/rc2 — 1-76. Deviations from the original BCS theory are 
understandable in that (i) the system is in the strong-cou pling regime, and (ii) the 
pairing is in a spin-triplet rather than a spin-singlet state. 1^ 

One measure of coupling strength is the ratio A(0)/|^| of the zero-temperature 
energy gap to the magnitude of the chemical potential. It is seen from Fig. [TTl that 
the strong-coupling regime is realized for / > 40 (i.e., A /i). At / = 20 the 
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system is in a transitional regime (A ~ fi). Another measure of coupling is the 
diluteness parameter n|ap, where a is the scattering length. In agreement with the 
first criterion, the matter is in the dilute (or strong-coupling) regime for / > 40, 
since this range corresponds to na^ — 0.63 < 1 when a is taken as the triplet n — p 
scattering length, 5.4 fm. A signature of the crossover from weak to strong coupling 
is the change of sign of the chemical potential, which occurs for / « 80 fFig. 
slightly below the crossover density between weak-coupling and strong-coupling 
regimes. 

In the limit of vanishing density, f —^ oo, the value of the chemical potential at 
T — tends to /i(oo) = — 1.1 MeV, which is half the binding energy of the deuteron 
in free space. Indeed, in this limit the gap equation reduces to the Schrodinger 
equation for a two-body bo und state, with the chemical potential assuming the 
role of the energy eigenvalue.E3 Thus, the BCS condensate of Cooper pairs in the 
'^Si-^Di state evolves into a Bose-Einstein condensate of deuterons as the system 
crosses over from the weak- to the strong-coupling regime. The crossover is smooth, 
taking place without change of symmetry of the many-body wave function. 

How does isospin asymmetry affect the transition? As the system is diluted, the 
critical asymmetry at which the pairing disappears increases from small values of 
the order of 0.1 up to the asymptotic value a — 1. The reason for this behavior is 
that in the low-density matter, the excess neutrons do not appreciably change the 
wave functions of protons, which are bound into pairs.l^ At asymptotically small 
densities, the chemical potential of protons tends to /ip(oo) = —2.2 MeV, which is 
just the binding energy of the system per half the number of particles bound into 
deuterons. The chemical potential of neutrons is determined by the excess particles 
in the continuum and goes asymptotically as ^„(cxd) (i.e. there is ultimately no 
energy cost in adding a neutron to the system). Note that the asymptotic behavior 
described is independent of the degree of isospin asymmetry. 

In closing this section, we call attention to the remarkable progress achieved 
during the last few years in trapping and manipulating ultracold fcrmion gases. 
The strength of the two-body interaction between the constituent fermionic atoms 
can be t uned using th e Feshbach resonance mechanism, by varying the external 
magnetic thus, the entire range from weak to strong couplings can 

be probed. Recent experiments on ultracold atomic gases have begun to explore 
their properties in cases where pairing occurs between atoms in different hyperfine 
states, which are unequally populated. ^'^ Systems of this kind are also subject to 
intensive theoretical study, wit h s pecific attention to homogeneous ASC phases^ 
phases with broken symmetrieJ^ and their realization in finite trap geometries'^. 
The universal features of ASC revealed by this effort should contribute significantly 
to our understanding of nucleonic pairing under isospin- asymmetric conditions. 
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5. Vortex states in compact stars 

5.1. Currents and quantized circulation 

The macroscopic physics of neutron-star rotation and its anomalies observed in the 
timing of pulsars can be described within the hydrodynamic theory of superfluids 
suitably extended to multifluid systems ^0*^ 103 104 105 rj.^^ elementary 

constituents of a Fermi superfluid - the Cooper pairs - are characterized by a co- 
herence length ^. On length scales L ^ ^, the condensate of Cooper pairs can be 
described by a single wave function tp, and the condensate forms a macroscopically 
coherent state. At an intermediate or "mesoscopic" scale, stellar rotation and the 
presence of a magnetic field lead to the formation of vortices, macroscopic quan- 
tum objects whose distinctive property is the quantization of circulation around a 
path encircling the vortex core. Since the condensate wave function must be single- 
valued at each point of the condensate, the circulation is quantized in units of 27rfi. 
On writing tp = ipQe^^, the gauge-invariant superfluid velocities can be expressed 
through the gradient of the phase of the superfluid order parameter x and the value 
of the vector potential A: 

Vr = -^VXr-—A. (46) 

In this expression, Cr = (e, 0) specifies the electric charge of protons (p) and neutrons 
(n) respectively, ttIt is their bare mass, and r stands for n or p. Applying the curl 
operator to Eq. 1)46(1 and implementing quantization of the circulation (with the 
phase of the superfluid order parameter changing by 27r around a closed path) , one 
finds 

curlt;,- = ^^Ur 5^'^\x — x^) = u^r , (47) 

j 

where ■nh/niT is the quantum of circulation, = ujr/^r is a unit vector along a 
given vortex line, Xrj defines the position of a vortex line in the plane orthogonal 
to the vector Ur, (5^^^ is a two-dimensional Dirac delta function in this plane, and 
B = curl A is the magnetic-field induction. The index j is summed over the sites of 
vortex lines. Eq. ((47|l treats the vortex cores as singularities in the plane orthogonal 
to Vr', this simplification is justified on scales larger than the coherence length of 
the condensate. For a single vortex, the integral of Eq. ((47|l completely determines 
the superfluid pattern. Since this equation is linear, the superfluid pattern created 
by a larger number of vortices is formed by superposition of the flows induced by 
each vortex. Obviously, the resulting net flow depends on the arrangement of the 
vortices. 

The condensate wave function can be written as ipi^) — /(^)e*^ in cylindri- 
cal polar coordinates (r,9,z). Upon integrating Eq. (|47|l . the neutron and proton 
superfluid velocities then become 
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where Ki is the Bessel function of imaginary argument. The divergence of the 
neutron- vortex velocity Vn{r) as r ^ is regularized by a cutoff A ~ The long- 
range nature of Vn{r) results in slow falloff of a density perturbation in the conden- 
sate. In a proton vortex, the supercurrent is screened exponentially on length scales 
of the order of the penetration depth A. Thus, for r :s> X, Ki{r/X) ~ exp(— r/A). 

On global, hydrodynamic scales, transition to a continuum vortex distribution 
can be carried through on the right-hand side of Eq. (|47f) by defining vortex densities 
72r = 'J2j {x — Xrj). Since the curl of Vn is simply 217 for rigid-body rotations, the 
number density of vortices in the neutron superfluid is related to the macroscopic 
angular velocity of the neutron condensate by the Feynman formula 

2m„r2 



(49) 



For typical pulsar periods P in the range 0.05 < P < 0.5 s, one has n„ ~ 
6.3. X 10^ P-^ ~ lOMO^ per cm^. In the case of a charged superfluid, Eq. (|47|l 
can be transformed to a contour integral over a path along which Vp ~ 0, since 
the supercurrent is screened beyond the magnetic field penetration depth A. If the 
proton superfluid is a type-II superconductor (i.e., A/^p > l/\/2), the continuum 
vortex limit leads to the estimate 

r) 

— ~ 5 X 10^^ cm-2 , (50) 

where $o = nhc/e is the flux quantum. We note that the number of proton vortices 
per neutron vortex is np/n„ ^ 10^'^ — 10^"', independently of their arrangement. 
The energy of a bundle of neutron or proton vortices is minimized by a triangular 
lattice with a unit cell area = (•\/3/2) d^. The lengths of "basis vectors" of the 
lattices in the neutron and proton condensates (the inter-vortex distances) are 



rf„ ^ ( — ] , dp^[^] , (51) 



where B is the mean magnetic-field induction. Using the estimates given in Eqs. (|49|l 
and H50|) . one finds that the neutron and proton inter- vortex distances are dn ~ 
10^^ — 10^'^ cm and dp ^ 10^^ cm, respectively. For typical values of the microscopic 
parameters, the penetration depth is of the order 100 fm = 10^^^ cm. Therefore the 
conditions ^„ ^ d„ and ^ min(A, dp) are satisfied, and the use of hydrodynamics 
on the local scale is valid. It is also clear that global hydrodynamics can be applied 
on scales that are much larger than dn (a fraction of millimeter). 

The strong interaction between the neutron and proton fluids gives rise to the 
entrainment efFcct: the supercurrent is a linear combination of the velocities of both 
fluidsf29 More specifically, the mass current and velocity vectors are related by 
a nondiagonal density matrix in isospin space. 



Pi \ / Pll Pl2 \ / Vi 
P2 J V P21 P22 J \V2 



(52) 
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where 1 and 2 label the isospin projections. The oflF-diagonal elements, which would 
vanish in the noninteracting limit, are evidently responsible for the entrainment 
effect. One fundamental conse q uence o f this effect is that the neutron vortex carries 
a non-quantized magnetic flusP^^^ of the same order of magnitude as the flux 
quantum $o. If the proton fluid forms a type-II superconductor, the number of 
proton vortices (sometimes called flux tubes) is 10^^ — 10^'^ per neutron vortex [see 
Eq. USUI) ]. Accordingly, the neutron- vortex motion (dynamics) is likely to be affected 
by the proton- vortex array. 

The arrangement of the proton- vortex lattice is a complex issue, and there exist 
several models for its configuration. 

(i) A class of (Hux-tube) models envisions the proton-vortex array to be spa- 
tially homogeneous, the vortic ity vector being inclined by some angle with 
respect to the spin vector.^^ Further, it is assumed that the flux tubes act 
as extended pinning centers for neutron vortices. In such models, a change 
in the neutron-vortex distribution is achieved by vortex creep of neutron 
vortices through the array of flux tubes. 

(iia) Vortex cluster models predict clustering of proton vortices over about 10% 
of the area occupied by a neutron vortex. One particular model generates 
a bundle of proton vortices coaxial with the neutron vorte x, th rough the 
entrainment currents induced by neutron- vortex circulation! I This gives 
rise to an average axisymmetrical magnetic field whose magnitude is com- 
patible with pulsar observations. 

(lib) The homogeneous distribution of proton vortices could be generically un- 
stable towards phase separation between a phase containing dense mesh 
of proton vortices and a phase devoid of vortices. A necessary condition is 
that the vortex lattice is sufficiently dilute, themean intervortex distance 
being much larger than the penetration depth.™ 

(iic) Proto- neutron stars are likely to possess natal magnetic fields; the nu- 
cleation of such a field will be associated with a first-order normal- 
superconducting phase transition, squeezing the field into bubbles of su- 
perconducting regions with high B ~ 10^^ G, and forming stable protonic 
vortex arrays which again cover about 10% of the total area! ^^*^ l The dy- 
namics of the neutron-vortex array in vortex-cluster models is controlled 
by the electromagnetic scattering of electrons off a vortex cluster. 

Current models of BCS pairing of protons do not exclude the possibility that 
there is a transit ion f rom type-II to type-I superconductivity of protons as the den- 
sity is increased !^ Type-I superconducting protons will have domain structures 
analogous to those observed in laboratory experiments on terrestrial superconduct- 
ing materials. The electrodynamics of the proton domain structures in NS can be 
treated by adapting the theories developed for laboratory superconductors, in which 
the magnetic fields are generated by normal currents driven around a cylindrical 
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cavity by temperature gradients.^EI A recent theoretical stud}^^ examined the 
effect of interactions between neutron and proton Cooper pairs on the status of 
proton superconductivity. The resuhs suggest that type-I superconductivity can be 
enforced throughout the entire stellar core if the strength of the interaction be- 
tween Cooper pairs is significant. However, within the mean-field BCS theory, the 
Cooper pairs are noninteracting enti ties, and any deviations from this picture must 
be due to fluctuations. Alford et al! I estimated the strength of the interaction 
between neutron and proton Cooper pairs due to fluctuations and found it to be 
too small to account for the interactions assumed in Ref. 11131 Nevertheless, type-I 
superconductivity of protons is not excluded by current calculations of ^Sq pairing 
of protons, and we shall address below its potential implications for the macroscopic 
manifestations of superfluidity in neutron stars. 

5.2. Constraints placed by neutron-star precession on the mutual 
friction between superfluid and normal-fluid components 

Since most of the inertia of a NS is carried by the neutron superfluids in the core 
and in the crust, the key to an understanding of NS rotational anomalies lies in the 
transfer of angular momentum from the superfluid to the normal (unpaired) compo- 
nent of the star, whose rotation is observed through the magnetospheric emission. 
At the local hydrodynamical scale, the rate of angular momentum transfer between 
the superfluid and normal components is determined by the equation of motion of 
a superfluid neutron vortex line. In the approximation that the incrtial mass of the 
vortex is neglected, the equation of motion is 

ujs{vs - vl) X 1/ + ({vl - vn) + C{vl ~ vn) X /y = , (53) 

where vs and are the superfluid and normal fluid velocities, vl is the velocity of 
the vortex, f is a unit vector along the vortex line, lus is the unit of circulation, and 
C, C are (dimensionless) friction coefficients, also known as the drag-to-lift ratios. 
These coefficients encode the essential information on the microscopic processes of 
interaction of vortices with the ambient unpaired fluid. Microscopic calculations 
commonly indicate ^' ~ 0, and one is left with a single parameter C. 

In the NS crust, neutron vortices are embedded in a lattice of neutron-rich nuclei, 
and the coefficient is determined by the interaction of the vortices with the nuclei 
and the electron plasma. (In some mod els, neutr on vortices are localized - pinned to 
the nuclei or situated in between them^E^Sllf the pinning is strong, Eq. H53|) is 
not valid, since the forces acting on the vortex are not linear functions of velocities. 
However, the regime of perfect pinning can be identified with the C — > cxd limit.) In 
the core of the star, the friction is controlled by the interaction of neutron vortices 
with the ambient electron-proton plasma; Eq. (|53() is valid under these conditions. 
Initial studies of the dynamical coupling between the superfluid and the normal fluid 
focused on interpretation of the observed post-glitch relaxation of pulsar rotational 
periods. It turns out that such interpretation is fraught with ambiguity, because the 
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long relaxation times can be obtained in the two opposite limits of weak (C — !■ 0) 
and strong (C — > oo) couplings. 

Recent observatiorP' of long-term periodic variabilities in PSR B1828-11, if at- 
tribut ed to precess ion of this pulsar, challenges existing theories of vortex dynamics 
in NS^^^^^^The importance of the inferred precession mode stems from the 
fact that it involves non-axisymmetric perturbations of the rotational state, re- 
moving the degeneracy with respect to ( that is inherent in the interpretation of 
post-glitch dynamics. In the frictionless limit, a star must precess at the classical fre- 
quency eil, where e is the eccentricity and is the rotation frequency. Clearly, then, 
there must exist a crossover from damped to free precession as C is decreased. The 
crossover is determined by the dimensionless parameters {Is/In)P and {Is/In)P' ■ 
Here, 13 = (/[{I - C')^ + C^], = 1 - /?(! - C')/C, Is is the moment of inertia of the 
superfluid, and In is the moment of inertia of the crust plus any component cou- 
pled to it oi l tini e scales much shorter than the precession time scale. The precession 
frequency i 



(54) 



In J In _ 

where fls is the spin frequency and e is the eccentricity. A no-go theorerrfrl^ states 
that Eulerian precession in a superfluid neutron star is impossible if {Is/In)C > 1 
(assuming as before C.' — > 0). There is a subtlety in this result: the precession is 
impossible because the precession mode, apart from being damped, is renormal- 
ized by the non-dissipative component of the superfluid/normal-fluid interaction 
((X /3'). In effect, the value of the precession eigenfrequency drops below the damp- 
ing frequency for any larger than the crossover value. This counterintuitive result 
cannot be obtained from arguments based solely on dissipation. In fact, according 
to Eq. H54(l . the damping time scale for precession increases linearly with C, and in 
the limit C — > oo one would (wrongly) predict undamped precession. If a neutron 
star contains multiple layers of superfluids, th e pic ture is more complex, but the 
generic features of the crossover are the same ! I While it is common to study 
perturbations from the state of uniform rotation, the precessional state may actu- 
ally correspond to the local energy minimum of an inclined rotor if there is a large 
enough magnetic stress on the star's core! ' ' 

In the core of a neutron star, the quantized neutron- vortex array is embedded in 
an electron-proton plasma, with the protons in a superconducting state. Electrons 
will scatter off the anomalous magnetic mon ients of (ungapped) neutron quasiparti- 
cles localized in the core of a neutron vortex.^2 Because of the spin-1 nature of 
the order parameter of the neutron superfluid in the core, the "^P^ vort ex core has an 
additional magnetization that scatters electrons more effectively.^^ An even more 
efficient scattering mechanism comes into play due to the flux ~ $o induced by the 
proton supercurrent on the neutron vortex via the entrainment effect If the pro- 
tons are non-superconducting in some regions of the core, then the strong nuclear 
interaction between protons and neutron quasiparticles localized within a vortex 
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core leads to a n efRcient coupling of the electron-proton plasma and the neutron 
superfluid.^^ The above models belong to the class of weak-coupling theories, i.e. 

<C 1, and, according to the no-go theorems, are compatible with free precession of 
the neutron star. However, these theories of mutual friction assume (unrealistically) 
that the proton- vortex lattice has no effect whatsoever on neutron- vortex dynamics 
in the core. 

The mechanism underlying mutual friction in the flux-tube models is the slow 
motion of neutron vortices through the pinning barriers (here, flux tubes) via ther- 
mally activated creep. Since, in general, the creep models presume that the vortex 
lattice closely follows the rotation of the pinning centers (N.B. the case of pe rfect 
pinning corresponds to C — > oo), the effective friction in these models is large ,1^^ 
C 3> 1. The kelvon- vortex coupling in the core provides another interaction chan- 
nel, leading again tdH^i C ^ 1- Similarly, for vortex-cluster models, in which the 
neutron-vortex lattice and the associated proton- vo rtex cluster move coherently, 
electron scattering by proton- vortex clusters also give 1. Accordingly, these 

theories are incompatible with free precession of a neutron star. This conclusion 
has been stressed by Link, -'^-'^^ and it has been argued that type-I superconductivity 
could be an alternative. 

In the crust of a neutron star, the neutron- vortex lattice is embedded in a lattice 
of nuclei and the charge- neutralizing background created by an almost homogeneous 
electron sea. In the vortex-creep models, the neutron-vortex lattice maintains ro- 
tational equilibrium via thermal and quantal creep through the pinning barriers 
(nuclei)! ' ■ ' I ^ 1 Hence these models imply C ^ 1 and are incompatible with free 
precession. If the pinning is a bsent, either because re-pinning cannot be achieved 
in post-jump equilibrium^^ and/or because of mutual cancellation of the forces 
from different pinning centers, the freely flowing neutron-vortex lattice interacts 
with t he elect ro n-phonon component of the crust. These interactions are weak and 
lead tcPaEHI 

C ^ 1. Note that the above estimates assume that the ratio Is /In 
is roughly of order unity. While it is difficult to estimate this ratio precisely, it is 
unlikely to differ from unity by many orders of magnitude. 



5.3. Type-I superconductivity in neutron stars 

The equilibrium structure of alternating superconducting and normal domains in 
a type-I superconductor is a complicated problem that depends on, among other 
things, the nucleation history of the superfluid phase. By flux conservation, the 
ratio of the sizes of the superfluid and normal domains is given by the relation 
ds/dN = \/IIcm/B ^ 10, where B ^ 10^^ G is the average value of the magnetic 
induction and Hem ^ 10^^ G is th e thermodynamic critical magnetic field. 
We first examine a 

modelini in which the magnetic field generated by the en- 
trainment effect supports the formation of domains coaxial with the neutron vortex. 
Consider a vortex that moves at a constant velocity Vl and carries a coaxial normal 
domain of protonic fiuid relative to the background electron liquid. Continuity of the 
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y 




Fig. 12. The structure of a rotational vortex placed in a type-I superconductor. The vortex 
velocity field is indicated by the concentric circles. The non-superconducting domain (shaded 
region) of radius a is coaxial with the vortex and carries a magnetic field Hem ~ 10^* G. The 
vortex motion along the x-axis generates a transverse electric field, which drives the electron 
current through the domain and causes Ohmic dissipation. 



electrochemical potentials of the superfluid and normal phases across the boundary 
between them entails the existence of a constant transverse electric field 

^^_m;v,v,{a) ^^^^ 
ea 

across the normal domain (see Fig. I12|l. The energy dissipated per unit length of 
a vortex is W = aE^ (a/b)'^, where a is the electrical conductivity and the factor 
(a/b)^ is the fractional area occupied by the domain. Combining this relation with 
Eq. H55|l . we obtain an alternative expression W = r/vj^ for the dissipation, which 



identifies the friction coefficient as 



129 



c = 



«ln('^^coth('^)-l 



A \ a 



(56) 



where $i — {P12/ Pii)^q- The zero-field conductivity of ultra-relativistic electrons 
is ctq = neC^cTc/peF, whcrc PeF is the electron Fermi momentum, Ug, is the electron 
number density, and Tc is the relaxation time for Coulomb scattering of electrons 
off protons in the normal domains. "'^^ The conductivity a = a^/ {lOcTc)'^ entering 
Eq. H56I) includes the effect of the bending of electron trajectories in the magnetic 
field. In this formula, ujc — eHcm/PF is the electron cyclotron frequency, which is 
proportional to the thermodynamic critical field Hem- For the typical density range 
p — (7.9 — 8.6) X 10^'' g/cm"^, the friction coefficient has the order of magnitude 
C ~ 10^^ at the temperature T — 10^ K and scales with temperature as C c>c T^. 
The drag-to-lift ratio satisfies the condition ( ^ 1 for all T below the critical 
temperature of the superfluid phase transition. 

We tu rn next to the friction coefficient as predicted by the model of Buckley 
et al.j^^in which the normal domains contain a large number ^ 10 neutron 
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vortices. In this case, the damping of the differential rotation between the electron- 
proton plasma and the neutron superfluid is due to the interaction of domain (i.e., 
non-superconducting) protons with the core quasiparticles confined in the neutron- 
vortex core. The relaxation process is therefore the same as in the case where the 
proton fluid is non-superconducting over the entire bulk of the core. However, the 
result needs to be rescaled by the ratio of the areas occupied by t he normal and 
superconducting layers. The relaxation time per single vortex i J I 



™ ^-^^exp 0.02^ , (57) 
\PFnJ nrriplanp \ epn-t J 

where ppp and ppn are the Fermi momenta of protons and neutrons, = 
m*m'^/(m* + m* ) is the reduced effective mass in terms of the proton and neu- 
tron effective masses, a„p is the total in-medium neutron-proton scattering cross 
section, A„ is the gap in the neutron quasiparticle spectrum, and epn is the neu- 
tron Fermi energy. [Eq. (|57|) differs from the analogous expression in Ref. 11241 by 
the factor 4m„//iP; here P is the pulsar period and to„ is the free-space neutron 
mass.] 

In the relaxation-time approximation and zero-temperature limit, the friction is 
written as 

PS^S PS^SCTnp 

where Up is the proton number density. For proton densities in the range pp = 
(4—8) X 10^^ g/cm^ and temperatures T ~ 10^~*, one finds C < 0.1. For a given 
model of the type-I superconducting structure, the friction coefficient C, must be 
rescaled by a factor [dN/dsY 0.01. 

Now we are in a position to discuss the implications of the friction mechanisms 
described above for neutron-star precession. The condition (Is/In)C < 1 seems 
to fulfilled unless Is /In ^ 1- The magnitude of the ratio Is /In depends on the 
superfluid/normal- fluid friction within all superfluid regions of the neutron star and 
is hard to assess. Glitches and post-glitch relaxation provide a model-independent 
lower bound. Is /In > 0.1. On the other hand, an upper bound is difficult to set. 
The deep interior of the star, if superfluid, could be decoupled from the observable 
parts of the star on evolutionary time scales, without any effect on short-timc-scale 
physics (although one does require C — > 0, rather than C oo, to prevent damping 
of the precession). At any rate, it is rather unlikely that Is /In exceeds unity by 
many orders of magnitude, and appealing to the lower bound on the ratio of the 
moments of inertia, one can conclude that the precession is undamped for both 
dissipation mechanisms considered. 

6. Concluding remarks 

This review has covered a number of aspects of nuclconic superfluidity, ranging from 
microscopic theories of pairing in nuclear systems and neutron stars to mesoscopic 
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frictional processes in superfluids and rotational anomalies in pulsars. Our survey 
of this important subfield is by no means complete. While the selection of topics is 
naturally biased toward the primary interests of the authors, we have chosen topics 
and problems with the intent of elucidating (i) the fascinating relationships between 
the physics involved at different scales and (ii) the richness of the contributions from 
diverse subfields of physics. We close by listing a number of issues and problems 
that call for further clarification and concerted effort within the general framework 
of our discussions. 

• The pairing problem at the level of mean-field BCS theory, with the pair- 
ing driven directly by in- vacuum nuclear interactions, is essentially solved 
within the density range over which these interactions are constrained by 
experiment. On the other hand, issues such as the screening of nuclear 
interactions, renormalization of the single-particle spectrum, and off-shell 
energy behavior of the pairing gap still defy quantitative resolution. Broadly 
speaking, extensions beyond the BCS theory are needed that incorporate 
fluctuation corrections while providing a consistent treatment of short-range 
correlations. 

• Superfluid phases with broken space-time symmetries have received much 
attention from theorists in recent years, while experimental realization of 

asymmetric superfluids opens the possibility of testing the predictions of 
theory. Importantly, relevant experiments are now probing the BCS-BEC 
crossover via the Feshbach resonance mechanism. Given a broad eff'ort, there 
is the prospect of mapping out the superfluid phase diagrams of interesting 
fermionic systems, both experimentally and theoretically, in the space of 
coupling strength, spin/isospin asymmetries, temperature, etc. 

• The rotational anomalies observed in neutron stars continue to provide use- 
ful constraints on the state of the superfluid matter in neutron-star interiors. 
Further theoretical studies of vortex dynamics, combined with pulsar timing 
observations, may be expected to shed new light on the internal structure 
of the superfluid phases of neutron stars, especially on the question whether 
protons form a type-I or type-II superconductor. 
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